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SUMMARY 


Methods Based on oblique- and normal-shock relationships and the 
continuity of mass flow through suitably chosen volume elements between 
the shock and body were developed to predict shock envelopes about two 
types of vehicles being considered for atmosphere entry. One type is a 
high-drag capsule shape. The other type is essentially a slender tri- 
angular wing capable of providing high lift or high drag, depending on the 
angle of attack. Predicted and measured shock envelopes were compared for 
a Mach number range of 3 to 15 for vehicles at high angles of attack; good 
agreement was found. Most of the available experimental data were in a 
speed and temperature range in which no Important real-gas effects 
occurred . 


INTRODUCTION 


During atmosphere entry, a vehicle is exposed to aerodynamic forces 
and heating rates which are related to the strength and disposition of 
the shock envelope about the vehicle. Present methods (e.g., refs. 1 
through 4) for defining detached shock envelopes are generally limited to 
simple shapes at zero or low angles of attack. Little information is 
available on predicting shock envelopes in the case of vehicles at high 
angles of attack. 

The purpose of this paper is to present methods for defining shock 
envelopes for the case of high angle of attack and to present a comparison 
of predicted results of the analysis with experimental results. Two types 
of vehicles, shown in figure 1, are considered. The vehicle on the left 
is a high -drag capsule shape, capable of providing low lift -drag ratios 
over a limited angle -of -attack range. The vehicle on the right is essen- 
tially a slender wing capable of providing low drag and high lift -drag 
ratios, and high drag and low lift -drag ratios, depending upon the angle 
of attack. These two types of vehicles will hereinafter be called capsule 
and conical types, respectively. 
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METHODS 


This section of the report will describe the general principles, 
relationships, and assumptions used in developing equations for predict- 
ing shock envelopes about capsule and conical types of vehicles; will 
present the derived equations; and will present the procedures and charts 
for applying the equations. Appendix A describes the symbols used. The 
detailed derivations of the equations are given in appendix B. 

Common to the analyses for both types of vehicles was the application 
of (l) the continuity of mass flow between the shock wave and the vehicle 
surface, (2) oblique- and normal-shock relationships, and (3) the unique 
correlation of stagnation -point velocity gradient with Mach number for 
different bodies which was pointed out in reference 1. The analysis 
follows that of Moeckel in the utilization of items (l) and ( 2 ). The 
analysis differs in that the sonic line is not involved except in its 
terminal position on the body and in that the shock is assumed circular 
in shape rather than hyperbolic . The circular shock is adopted in the 
interest of generality since it is a simple compromise between the 
n hyperbo lie -like” shock about cylinders or spheres and the Tl elliptic -like" 
shocks about long plates or discs. Moreover the shock dimension is then 
characterized by a single value R thus simplifying the analysis. 
Application of item ( 3 ) gives a means of determining the value for R as 
will be shown. The gas properties required in these applications were 
taken to be those of perfect gases or real gases in equilibrium. The 
analyses for the two types of vehicles differ as to detail and thus are 
described separately in the following two sections . 


Capsule -Type Vehicles 


Figure 2 illustrates the mass -flow components that are taken into 
account in the analysis of capsule-type vehicles . It is assumed that 
the shock wave is composed of circular-arc elements whose radii and loca- 
tion with respect to the body are to be determined. The sonic -point 
locations on the body are based on experimental and theoretical results 
(e.g., ref. 2 ) which indicate that for spheres, sonic points exist at 
surface elements inclined at about ^5° to the stagnation streamline . In 
cases for which the spherical segment ends before the surface inclination 
attains V? 0 , the sonic point exists at the edges of the body. 

The foregoing considerations and the application of mass -flow 
continuity resulted in the following equations: 
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Equation (l) relates the ratio of the sonic-point shock standoff distance 
to shock radius, A/R, to the shock angular position 0 opposite the sonic 
point, the density ratio across the shock at that position, Pq / P 0 > and the 
ratio of the mass flow per unit area at the sonic point, p*V*, to entering 
mass flow per unit area P o V 0 . Equation (2) is a geometric identity 
relating the sonic standoff distance to the standoff distance at the stag- 
nation point, the hody radius, r, and angular position of the sonic point, 
9, and the shock radius and angular position, 0, opposite the sonic point. 
This equation is strictly applicable only to the axially symmetric case. 

In the general case, involving angle of attack, the stagnation streamline 
is curved and the expression relating A and A 0 is more complex than 
equation (2). For many cases, however, the curvature is small so that 
» length A 0 is substantially the same as its horizontal projection. The 
approximation is made, therefore, that equation (2) is applicable to the 
general case . 

The foregoing equations based on continuity and geometric relation- 
ships are not sufficient for a solution of shock -wave shape and location, 
because of the existence of three unknowns, A, R, and 0, and only two 
equations. The remaining quantifies are either known (r, cp) , calculable 
independently of the unknowns ( P*V^/poV 0 ) , or calculable as functions of 
the unknowns (A 0 , Pg/ P Q ) • An independent equation giving additional 
relationships among the shock parameters is therefore needed. The third 
equation was derived on the basis of the previously mentioned velocity 
gradient correlation. The equation that follows defines a function which 
correlates the standoff distance at the stagnation point, the body stagna- 
tion point and shock radii, and the density ratio for various axisymmetric 
shapes : 


f 




(3) 


The values of the correlation function f as a function of density 
ratio are shown in figure 3- The recommended values indicated were 
obtained from the derived correlation formula by using theoretical stag- 
nation point standoff distances and shock radii given in reference 2 for 
spheres . 


1 + 


To define the shock radius and shock location, the three equations 
are solved simultaneously. For a solution it is necessary to relate 
p 0 / p 0 with Pl /p 0 and 0 in equation (l). This relationship depends on 
the nature of the gas involved. For a perfect gas, Pg/p 0 can he expressed 
as a function of 0 and p 1 /p Q . The latter ratio and the other flow 
parameters, p*V* / p 0 V 0 and P s /P 0 , can "be: calculated from perfect -gas 
relations . In the case of a real gas , charts or tables of the thermody- 
namic properties and normal -shock properties are required (refs. 5 and 6, 
respectively, in the case of air). 


As a convenience in obtaining solutions of shock envelopes for the 
case of real air, the charts of figure 4 were prepared. These charts 
present values of A/r and 0 as functions of the normal -shock density 
ratio, Pj/pq, and cp. The application of these charts in determining 
the shock envelope is as follows: Consider the case of a vehicle at angle 

of attack, as shown in figure 2. We will confine our attention to the 
trace of the shock envelope on the vehicle's vertical plane of symmetry. 
These segments have radii R x and R 2 which coincide tangentially at the 
origin of the coordinate system X-Y. To locate the body with respect to 
these shock segments, the coordinates of the two sonic points (l) and (2) 
as well as the shock radii must be found. These values are given by the 
following equations: 
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The values of 0 and a/r are located from the charts of figure 4 at the 
appropriate values of shock density ratio Pi/p Q an( i sonic -point incli- 
nation 9 . The values (X 2 - Xx) and |Y 2 | + i Y x l are the known distances 
between the sonic-point coordinates (l) and (2). An example calculation 
is given in appendix B. 
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In preparing the charts presented in figure 4 , the following three 
approximations were made: (l) for a value of Pi/ P Q greater than 5, the 

value of Pg/p 0 was assumed to he equal to Pi/ P 0 j ( 2 ) the value of 
P*V*/ Po V 0 was always assumed to be equal to that given by the perfect - 
gas relations with 7 = l.b (eq. (Bl)); and (3) the value of p s /p 0 was 
always assumed to be that given by the perfect -gas relations with 7 = 1.4 
(eq. (B 7 ))- The possible errors introduced by these approximations were 
assessed in the following manner. 

Real-gas effects on p s / p Q and P*V^/ p o V 0 were assessed by comparing 
the ideal-gas values (see eqs . (Bl) and (B 7 )) with the values for equi- 
librium flow as calculated with the aid of the charts of reference 6 at 
various levels of density ratio, pressure, and enthalpy. It was found 
that no large differences existed between the values representing ideal ■ 
and real gases . Small differences can be expected because of the weak 
dependence of p^V*/ PpV Q and p Q / p s on the specific heat ratio 7. In 
view of the foregoing, significant errors in the calculated shock dimen- 
sions then can arise only through the effects of large departures of the 
actual ratio Pq/Pq from idealized conditions. To assess such errors, 
this effect of a real gas on the shock parameters 0, A, and A 0 was 
evaluated. Values of these parameters for a sphere in air in equilibrium 
flow at a velocity and absolute density at which a large variation in 
p q/ p 0 occurs were compared with values obtained in a hypothetical gas 
having the same normal shock density ratio but with Pq/ P Q independent 
of 0. The calculated differences in A 0 were negligible, but a small 
effect on shock radius resulted in a difference of about 15 percent in 
A* between air in equilibrium and the hypothetical gas. 


Conical-Type Vehicles 


Figure 5 illustrates the mass -flow components that are taken into 
account for vehicles of conical type. Free-stream mass -flow components 
exist both in a crossflow direction (p 0 V 0 sin a 1 ) and in a tangential 
direction ( p 0 V 0 cos a f ) as indicated. The shock-wave trace in the cross - 
flow plane is taken as an arc of a circle. Sonic locations on the body 
are assumed to exist at points where the body surface inclination with 
the crossflow plane is ^ 5 °, Tor all finite ratios of a/b, and at the edge 
for the case of a flat plate (a/b = 0 ). For an infinitely slender ellip- 
tic cone, the determination of the shock envelope is essentially a two- 
dimensional (crossflow) problem. The shock dimensions in such a case are 
determined by the following equations which are analogous to equations 

(1) to (3). 
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A 
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a/b 

J 1 + (a/b) 2 _ 


- (l - cos 0) 
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Pq 2 + (Ap/ r) 
Ps _1 + (Aq/ r) 



( 6 ) 

( 7 ) 


In cases where the cone is not infinitely slender, the tangential 
flow component must also he taken into account. Experimental evidence 
obtained for the present study indicates that the shock envelope is con- 
ical up to high angles of attack. With the shape of the shock thus 
defined, it is possible to derive expressions on the basis of flow conti- 
nuity and the geometric properties of the body to determine the shock 
standoff angle P as well as the shock radius R indicated in figure 5* 
The details of this derivation are developed in appendix B and give the 
following : 
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a ! = 3 + a + e z 

Values of f used were the same as those for the capsule. The basis for 
this is explained in appendix B. Values for 0 and A/r are given in the 
charts of figure 6 as a function of density ratio Pj_/ P Q for various 
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parameters a/b . The same assumptions were made in regard to gas proper- 
ties as in the charts of figure k. In the application of Q and a/r in 
equation (8), however, the density ratio is PiVpq* that of an oblique 
shock at angle a* . Since a* is a function of 0, the angle in question, 
iterative calculations beginning with an initial assumption for the value 
of 0 are required. However, trial computations quickly converge to the 
solution. The validity of the solution is - restricted at high angles of 
attack to the angle range wherein the sum of the thickness angle e z , the 
shock angle p, and the angle of attack a is less than approximately 
90°. When the sum of these angles is greater than 90°, the tangential 
flow component is reversed and the flow can no longer be conical. The 
validity is restricted at low angles of attack; if Mc>sin a 1 is less 
than about 2, the components of mass flow are not given accurately in the 
present simplified treatment. 


COMPARISONS OF ANALYTICAL AND EXPERIMENTAL RESULTS 


This section of the report will first establish the validity of the 
results of the present method in determining shock standoff distance for 
simple symmetrical shapes at zero angle of attack by comparison with the 
results of the more accurate methods of Van Dyke and Belotserkovski and 
with those of experiment. Next, the results of the method in determining 
shock shapes of capsule- and conical-type vehicles at angles of attack 
will be compared with experimental measurements . 

Figures 7 and 8 show a comparison of analytical results for symmet- 
rical shapes at zero angle of attack with available experimental results . 
Agreement between the results of the present analysis with those of Van 
Dyke and Belotserkovski and with experimental value's is satisfactory. 

Figure 9 gives examples of experimental and estimated shock envelope 
traces about capsule -type vehicles for several angles of attack and Mach 
numbers. Only one example is presented wherein appreciable real -gas 
effects were involved. This is the capsule at Mach number 1^.9 which was 
tested at a simulated altitude of 30,000 feet at a velocity of 10,l8o 
feet per second. Generally, the experimental and estimated shock envelopes 
are in good agreement . 

In figure 10 calculated and experimental results for the shock angle 
of a family of elliptic cones are shown as a function of angle of attack. 
The cones have thickness -to -span ratios of 0 (a flat plate), l/3> and 1 
(a circular cone). The results are for a Mach number of b. The limiting 
lines shown indicate the angle -of -attack range for which the analysis is 
applicable. Agreement between predicted and experimental results is good. 


In figure 11 predicted and experimental shock standoff angles are 
compared for "blunted as well as sharp-nosed conical shapes . The nose of 
a blunted conical shape is treated as a capsule-type vehicle and the 
shock angle for the conical part of the shock envelope is faired to the 
capsule solution at the point of equal slope. The data of figure 11 are 
for angles of attack from 20° to 70° and for Mach numbers from 3 to 10. 
Good agreement between experimental and theoretical results is noted. 


CONCLUDING REMARKS 


Simple methods were developed for predicting shock envelopes about 
two types of atmosphere entry vehicles at angle of attack. Predicted 
shock envelope shapes for capsule -type vehicles were compared with exper- 
imental values for air flows within the Mach number range 3 to 1^.9 and 
the angle -of -attack range up to 33 ° • Predicted shock standoff angles for 
conical vehicles were compared with experimental values for air flows 
within the Mach number range 3 to 10 and the angle -of -attack range 20° to 
70 ° . 


The methods are applicable to perfect -gas flows and to equilibrium 
flows of real gases. Charts were developed specifically for predicting 
shock envelopes . 

Satisfactory agreement between predicted and experimental shock 
envelopes was found for the range of flow conditions covered by the com- 
parisons . Except for one case, only perfect-air conditions were covered, 
thus essentially limiting the evaluation of the methods of prediction to 
the perfect -air case. For the evaluation of the methods for high -enthalpy, 
equilibrium flow of real gases, further experimental research is needed. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Mar. 6 , 1961 
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APPENDIX A 
NOTATION 

a vertical axis of elliptic cone, ft 

b horizontal axis of elliptic cone, ft 

b* spanwise distance to sonic point on elliptic cone, ft 

f shock correlation function 

M Mach number 

m mass -flow component 

R shock radius , ft 

r stagnation point body radius, ft 

V stream velocity, ft/sec 

X,Y,Z Cartesian coordinates 

A shock standoff distance from sonic point on body, ft 

A 0 shock standoff distance from stagnation point on body, ft 

a angle of attack, deg 

a’ conical shock angle with respect to free-stream direction 

in vertical plane of symmetry, deg 

3 angle between conical shock and body in vertical plane of symmetry, 

deg 

5 flow deflection angle through oblique shock, deg 

e half -angle subtended by capsule forebody arc, deg 

€y cone semiapex angle in plane of horizontal axis, deg 

e z cone semiapex angle in plane of vertical axis, deg 

0 inclination of shock element opposite sonic point on body with 

respect to free stream or to crossflow stream (see fig. 2), deg 



specif ic -heat ratio 
density, lb/ft 3 

density behind normal shock, lb/ft 3 

density behind shock inclined at angle a 1 , lb/ft 3 

inclination of body surface at sonic point with respect to plane 
normal to free -stream direction, deg 

Subscripts 


free -stream conditions 
stagnation point on body 

point behind shock element at 0 degrees 


sonic point on body 
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APPENDIX B 

DERIVATION OP EQUATIONS 


The following sections give the details of the derivation of the 
equations defining the shock envelopes for capsule- and conical -type 
vehicles. Numerical examples are also presented. 


CAPSULE -TYPE VEHICLES 
Flow. Continuity Equations 


Mass -flow continuity is satisfied for an elemental volume between the 
shock and the body as shown in figure 12. The flow is three-dimensional 
and assumed to proceed radially from the stagnation streamline . Mass 
flow passes out from the volume through a strip of unit width and length 
A. The arithmetic mean of the vertical components of mass flow at the 
sonic point on the body and at the shock -wave position, 0, opposite the 
sonic point is utilized. The component at the sonic point is equal to 
the product of cosine cp and the mass flow passing through a normal shock 
and expanding isentropically to sonic value. This product, expressed as 
a mass flow per unit area and normalized with respect to the free -stream 
value p n V n , is _ _ x 

7-1 

cos cp (Bl) 

The component of mass flow per unit area at the shock normalized with 
respect to the free -stream value is given by oblique -shock relationships 
as 



Hl- 1)sln 20 


The mass flow passing out of the volume element is then 

a r p o v o /pa \ 

f ^ (a ■ V sin 20 + p * v * cos 9 


(B2) 


(B3) 
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Similarly, the product of entering unit -area mass flow p Q V 0 and the 
projected "wedge” strip area (l/2)R sin 6 gives 

S2X2.R sin 9 (B4) 

2 

Expressions (B3) and (b 4-) are equated and rearranged to give equation (l). 


A 

R 



( 1 ) 


Geometric Relationships 


A necessary step in. subsequent shock solutions is to relate the shock 
standoff distance, A, at the sonic point to the standoff distance, A 0 , at 
the stagnation point. Since the capsule forehody and the shock elements 
are circular, these standoff distances are related as follows: 

A = A 0 + r(l - cos cp) - R(l - cos 0) (B5) 

The above equation normalized with respect to shock radius R gives 
equation (2). 

— = — + ~ (l - COS 9) - (l - COS 0) (2) 

R R R 


Velocity Gradient Correlation 


is a function of Mach 


It has been shown in reference 1 that the stagnation -point velocity 

gradient in the dimensionless form — — — - — is a function of Mach 

V 0 (1 + A 0 /r ) 

number (or normal -shock density ratio) and correlates data involving a 
wide variety of axisymmetric shapes . In the following derivation the 
above stagnation -point velocity gradient will be equated to a more conven- 
ient correlation function in which the velocity terms are eliminated and 
in which only shock standoff distance, body radius, and shock radius terms 
appear . 


Continuity of mass flow is considered for a small conical volume 
element centered on the stagnation streamline as shown in figure 13 • The 
line 0-1 represents an element of the volume extending from the shock to 
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the "body. The flow is deflected (5 + dcp) with respect to 0-1 at point 0 
on the shock. Oblique shock relationships for small 6 give 
5 = (pj/Pq - l)d0. Geometrical considerations give (A 0 + r)dcp = R d0, so 

\ A 0 + r 


that & + dcp 
normal to 0-1 is then 


= [7£i.- 

LVPo 


/ 


PoVc 


+ 1 




dcp* The unit -area mass flow at 0 


+ 1 


dcp. The unit -area mass 


flow at 1 is [d(pV)/ds]r dcp, a value to be determined. The arithmetic 
mean of the mass flow leaving the conical surface of length A 0 is equated 
to that entering the projected surface of radius (r + Ao)dcp. The resulting 
equation is 


Y | 2lt (r + A 0 )dcp(p 0 V 0 ) ^ 


A 0 + r 


+ 1 


dcp + 2jtr dcp 


' a(pv) 

ds 


r dcp 


= P 0 V 0 «( r + Ao) 2 dcp2 


The stagnation -point mass -flow gradient d(pV)/ds given by the above 
identity is 


a(pv) _ Pq v o (- 1 + &o\ 
ds A 0 \ r J 


1 





(B6) 


At the stagnation point V = 0, p = p g> dp/ds = 0, so d(pV) = p s dV. 
Equation (B6) is rearranged to give 


where 



(B7) 


It is convenient, for reasons to be discussed later, to apply the factor 
2 to equation (B7), thus 



or, letting f replace the left-hand term, 


f = 2 £ 2 . 1 - (1 +£2 
p s L \ r. 


-Ot 


Numerical Example 


The charts of figure 4 are now applied to an example calculation of 

the shock trace about the lenticular forebody shown in figure 2. The fore 

body arc subtends a half -angle e of 30°. The flight condition is 
M = 6 . 3 * Pi/ p 0 = 5*33 j angle of attack of 1 6 °. At this angle the upper 

sonic point (l) is at the corner and is inclined cp x = ( 30 ° - 1 6°) = l4° . 

The lower sonic point (2) is on the body where cp 2 = 45°. The relative 

locations of the sonic points are (x 2 - X x ) = 4,30 and lYj + |Y 2 | = 15.60 
The values A/r x = 0.106, 0 X = 16*5°, A/r 2 = 0.175> and 9 2 = 31.4° are 
found in figure 5* Equations (4) are numerically evaluated as follows: 


15.60 - 4.30 ( — 2^521 — \ 
Xi = VO. 175 + 0.14 77 

f 0.284 3 + 7 0.521 

Vo. 106 + o.o 4 iy Vo. 175 + 0.147: 


v _ 2.43 x 0.284 _ ), ^ 

i x ~ 'i ' * f 

0.106 + o.o 4 i 


Yp = 1 


Ri = - 


g Ha 
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CONICAL VEHICLES 

Flow Continuity and Geometric Considerations 


Detailed derivations of equations defining the shock envelope about 
conical vehicles are developed with the aid of figure l4. In this figure 
are shown three views of a section of an elliptic cone at angle of attack 
a. The flow is assumed to be conical. The shock element in the vertical 
plane of symmetry inclined at angle a r is taken to be along the 
X axis . A volume element of length dx between the shock and the body 
surface is bounded by the vertical plane of symmetry poo r p* and the 
surface containing sonic line nn’ and its projection onto the shock 
surface, mm* . Continuity of mass flow is satisfied for this volume ele- 
ment . 


To write the continuity equation it is necessary to evaluate projected 
areas of the various surfaces of the volume element . Area mnop is 
similar to area m’n’o'p’ since the shock envelope and body surface are 
conical with respect to a common origin. Thus 


or 


where 



^m’n’o’p’ 
(I + dxf 


Sm*n t o , p t 


Snrnop 





> (B8) 


Smnop 


y (2A 0 + A) 


It should be noted that the control surface intercepting the tangential 
mass -flow component is not normal to the axis of the body. Therefore , 
in the plane of the control surface* the elliptic ity of the associated 
body section is different from the true cross-sectional value. The 
relative distortion in ellipticity is of the order cos(e z - P). Since 
€ z and (3 are usually small and in any event tend to cancel* the error 
resulting from section distortion is neglected. The distance from the 
sonic point (where surface slope of ^5° is assumed) to the vertical plane 
of symmetry is 



R sin 0 = b* 


(B9) 


Equations (b 8) to (BIO) and the terms indicated on figure 1^, are used 
to project the following areas on the Y-Z plane. 


b(3A + 2Nb) 
Smnop , == 


3 A + If 


R(l - cos e)h'dx b 2 (sin 0)dx 

* mrvn * — ’ * ' ' - - - - - r -r ■ ' 1 1 1 

3* 


Sm^pp 


1 + l £ 

sP 


Sm’n* 


Ab dx 


nm 


1 + ( 1 l 


on the X-Y plane , 


Snm ! p r p 


b r dx 


b dx 



and on the X-Z plane , 


Smnn *m 1 ~ A dx 

Next, the values of the mass -flow components per unit area ( pV ) 
normal to each area are required. All velocities in the X direction 
are taken as V Q cos a T , being substantially tangent to the shock. The 
density within the control volume varies from a maximum of p s ! > Pi 1 
at line oo* to a minimum of p* < p x T at the sonic line n n f . The 
density distribution is approximately accounted for by assuming p f ai 
the mean density of the flow through the surfaces mnop and m , n , o t p t , < 
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p* as the mean density through the surface m’n’nm. The sonic value 

p* = Pi’ [2/(7 + l) - (7 - l)( P 0 /Pi‘ is used. The unit-area mass 
flow (m) in the Y direction through mnn'm* is p 0 V 0 [sin q/(a/r)]. 
Values for Q and A/r are determined from figure 7 with Pi’/p 0 in place 
of Px/pQ' Summation of all unit mass -flow products with their associated 
areas is then made to satisfy the continuity of flow. Omitting the 
detailed algebraic steps involved yields the following equation: 


A 

b 


tan €y 

sin e _ 4 Pi* N 


tan a’ 

.6 *il + (a/b) 2 3 P ° J 




+ tan e y (2 - BL.) 

tan a* \ Pi/ 



(Bll) 


Equation (Bll) along with the following geometric relationships 
permits calculation of the shock angle j3 as a function of angle of 
attack a . 


0 


tan ” 1 



tan e-y 


(8) 


and 



a 1 = p + a 4 - e z 


Velocity Gradient Correlation 


Equation (7) is derived on the basis of two-dimensional flow through 
a differential volume element in the vicinity of the stagnation streamline 
in a manner similar to that in the case of axisyrametric flow. Values for 
f given by equation (7), plotted in figure 3 > were obtained by using 
the theoretical stag nation -point standoff distance and shock radii given 
in reference 3 for circular cylinders . Under conditions of the same 
ratio for A 0 /v o the stagnation-point velocity gradient for axisymmetric 
flow can be shown to be one -half the value for two-dimensional flow due 
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to radial M relief. n Accordingly, the factor 2 applied to the axisymmetric 
velocity gradient places it in correlation with two-dimensional gradients, 
thus relating both two- and three-dimensional shock distance parameters 
to the same function f of density ratio. 

The derivation of two-dimensional stagnation-point velocity gradients 
proceeds as follows: The line 0-1 in figure 13 in two-dimensional flow 

represents the side surfaces of a wedge of unit depth. The mass flows 
per unit area at 0 and at 1 remain the same as in the axisymmetric case. 
They are, respectively, 


Po^o 




Aq + r 
R 


+ 1 


dcp 


and 


iiellr a? 


The product of the entering mass -flow per unit area with the area of 
surface 0-0 is p 0 V 0 2(r + A 0 )dcp. This is equated to the product of the 
mean unit mass flow leaving surface 0-1 with the total area 2A 0 . 





Aq + r 
R 


+ 1 



= P 0 V 0 2( r + A 0 )dcp 


The stagnation -point mass -flow gradient d( pV)/ds is given by reducing 
the above identity to 


d( pV ) _ PqVq /g + Aq\ 
ds _ A 0 |_\ *V 


G + t ) (k - 0 A t 


(B12) 


At the stagnation point Y = 0, p = p s , dp/ds = 0, so d(pV) = p s dV. 
Equation (B12) is arranged in the form 


or 



f 







( 7 ) 
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It is to be remembered that r is the stagnation -point body radius; 
therefore, the application of equation (7) for elliptic bodies requires 
b 2 /a in place of r. 


Numerical Example 


The shock angle (3 is calculated for an elliptic cone with €j = 15°, 
e z = 5°j and a/b = l/3 at angle of attack a = ^5° at flight Mach number 6. 
An initial assumption that (3 = 5° is made; thus a trial value for 
ad = 5° + ^5° + 5° (see eqs . (8)). The crossflow Mach number is 
6 sin 55° = ^*91 a t which p 1 l /p Q = ^*97- Corresponding to this density 
ratio with a/b = l/3, the charts of figure 6 give A/r = 0.23^ and 
0 = I3.6 0 . Equations (8) are then numerically evaluated as follows: 


N 





1 - 0.972 1 

3--M0 8 


0.237 + 3 



J 




-0.117 


A 

b 


i , 0 . 268 

P • it ( 4 . 97 )( - 0 . 117 ) 

1I27 

6 jl + (1/3 f 

J 1+ (s) 

70 . 237 ) + 0.268 ( . 0.693,^.97 

Vo. 23 V 1A27 


0.50^ 


(3 = tan- 1 [(0.W - 0.117)0.268] = 5. 9° 


The calculated value for (3 is sufficiently close to the assumed value 
that an iterated calculation is unnecessary. 

The above example utilized the result for Qj 1 / 9 0 as "by oblique 

shock relationships for air with 7 = lA. At higher velocities and at 
altitudes where large departures from ideal gas properties occur recourse 
must be taken to shock density versus velocity and altitude plots such as 
are given in reference 5* -An initial value for PxVpo determined 
from such plots as that corresponding to the velocity V Q sin ad at the 
altitude in question. Generally, one iteration is sufficient to fix ad 
and thus p 1 t /p Q . 
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Figure 3*- Shock correlation function. 
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Figure 4 -.- Charts of A/r and 9 for capsule -type vehicles. 
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Normal shock density ratio, 

Figure 6.- Charts of A/r and Q for conical -type vehicles. 






Figure 7-“ Shock standoff distance for two-di 
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Figure 11.- Shock angles for conical i 
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